Non-Markovian effects in the quantum noise of interacting nanostructures 
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We present a theory of finite- frequency noise in non-equilibrium conductors. It is shown that Non- 
Markovian correlations are essential to describe the physics of quantum noise. In particular, we show 
the importance of a correct treatment of the initial system-bath correlations, and how these can be 
calculated using the formalism of quantum master equations. Our method is particularly important 
in interacting systems, and when the measured frequencies are larger that the temperature and 
applied voltage. In this regime, quantum-noise steps are expected in the power spectrum due to 
vacuum fluctuations. This is illustrated in the current noise spectrum of single resonant level model 
and of a double quantum dot -charge qubit- attached to electronic reservoirs. Furthermore, the 
method allows for the calculation of the single-time counting statistics in quantum dots, measured 
in recent experiments. 
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INTRODUCTION 

Vacuum fluctuations are one of the most intrigu- 
ing consequences of the quantum theory. In elec- 
tronic systems, they manifest as electron- hole cre- 
ation/annihilation processes in a time given by the 
Heisenberg uncertainty relation, t ~ 1/cj, being uj the 
measuring frequency. In order for these processes to be 
seen, other types of fluctuations must be overcome. For 
example, a system in thermodynamic equilibrium must 
be at a temperature T much smaller than this frequency, 
and in a system driven out of equilibrium, such as a 
mesoscopic conductor subject to an applied voltage V, 
the quantum- noise regime (QNR) reads Hlj ^> ksT^eV. 
Zero-point fluctuations in quantum-transport systems 
were first measured by Schoelkopf and collaborators pQ 
through the current-noise spectrum [2] 

/oo 
dre- <WT <{J(t + r),J(t)}) c , (1) 
-OO 

which reveals valuable information beyond that con- 
tained in the dc current [3HS]. Among the various 
methods to calculate S^ 2 \oj) : quantum master equations 
(QMEs) are particularly attractive because of their sim- 
plicity and generality for treating dissipative dynamics 
of interacting systems [5H2]- Typically, the Markovian 
approximation (MA) in the system-reservoir coupling is 
employed. This, however, fails in describing the noise 
spectrum in the QNR [10 , and although there have been 
a few attempts to go beyond the MA in the context of 
QMEs [11 , a complete noise theory is yet lacking. 

In this paper we present such a theory. Our method al- 
lows the calculation of the current and voltage noise spec- 
trum of a system described by a generic non-Markovian 
QME, and can be applied to the increasing number of ex- 
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FIG. 1: (color online). Schematics of counting: The density 
operator evolves from the initial separable state at time to 
(represented by two distinct ellipses) until it reaches a steady 
state at time t = 0, where it is no longer in a product state 
(single ellipse). At time t = counting begins. The shading 
highlights the time interval where counting is effective. Full 
circles denote tunnel vertices with counting factors x 0> 
empty circles denote standard tunneling vertices (% = 0). 
Contractions between tunneling events in counting and non- 
counting intervals (dashed over-line) give rise to T(x, z), while 
contractions within the counting interval (solid over-line) give 
rise to the self-energy H(x, z). 

periments exploring the QNR [12 j. The theory naturally 
contains the physics of vacuum fluctuations, for which a 
proper inclusion of initial system-bath correlations is es- 
sential. Furthermore, the method enables to determine 
the charge-noise spectrum 

/DO 
dre-^({Q(T),Q(0)}) C) (2) 
-OO 

as it is shown for a single resonant level (SRL) model. 
This noise dictates the back-action when the conductor 
is used as a detector of another quantum system [T3] . 
The technique is used to study the full noise spectrum 
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of a double quantum dot charge qubit in the hitherto 
unexplored QNR. As we will see, in this regime transport 
fluctuations are mediated by the zero-point dynamics, 
showing a series of steps at frequencies corresponding to 
resonant processes in the system. 

THEORY 

Here we consider phenomena that can be described by 
the general QME 

p(t) = Cp(t), (3) 

where C is the Liouvillian, that governs the evolution 
of the density operator (DO), p, describing the dynam- 
ics of the total system. Specifically, we focus on the 
case in which a central system exchanges particles with 
a bath, and this exchange is amenable to the count- 
ing of particles. We will take here the case of trans- 
port through a central quantum coherent system, at- 
tached to fermionic contacts. The Hamiltonian of the 
system is of the form H = %s + + Hy. Here 
Hs = ^ a Ea\ a )( a \ is the central-system Hamiltonian, 
with E a the energy of the 7V a -electron many-body eigen- 
state | a). The left/right reservoirs (at equilibrium with 
chemical potentials Pl/r — E F ± eV/2) are described 

by u R = E fe . a£L,R £ kctCka c ka: with Ska the energy 
of the k-th mode in lead a. The tunnelling Hamil- 
tonian is given by U v = Y^kam v kamcl a d m + H.c. = 
Efea m Ea )a ' V kam cl a (a\d m \a')\a)(a'\ + H.c, where c[ a 
creates an electron with momentum k in reservoir a and 
d m is the annihilation operator for the single-particle 
level m in the central-system. Vkam is a tunnelling am- 
plitude and e = h = ks = I throughout the text. Un- 
der the previous Hamiltonian, the DO evolves accord- 
ing to equation (J3|, with £• = — i [Hs + T~Lr + Hv, •] = 
(Cs + £r + £v) # - We are interested in the central- 
system dynamics, for which we consider the reduced sys- 
tem DO ps(t) = TVr {/)(£)}. If we choose to to be the 
time at which system and reservoirs are in a separable 
state, p(to) = ps(to) ® Pr\ with ps(to) arbitrary and p^ 
the equilibrium bath state, the evolution of ps(t) in the 
Laplace space is given by 

Ps(z) = TrR | [z - C}- 1 p s (t ) ® p e R q } = n (*)ps(*o)- (4) 

Here, we find the propagator Qq(z) = [ z — W(z)] -1 , with 
kernel W(z) = Cs + Y,(z), being Y,(z) the non-Markovian 
(NM) self-energy, and whose form can be derived using 
the expansion 

- J T = / r jtUv / r V. (5) 
z-C z-Cs-C R ^ o \ z-Cs-C R J 

This gives 

S(z) = Tr R ley - _ 1 _ £ £vf%] + ■■■ (6) 



Technical details on how to evaluate this expression [TH - 
US] are not relevant for the main discussions and are given 
in appendix]] 

Cumulant generating function 

Our goal here is, given Eq. Q, to derive a formula 
for the cumulant generating function (CGF) in terms of 
known quantities such as the self-energy. This will al- 
low us to calculate NM current correlations up to arbi- 
trary order at zero frequency. Furthermore, we aim to 
give an expression for the NM finite-frequency noise cor- 
relation function. If the transfer of electrons between 
system and reservoirs is amenable to counting, the full 
counting statistics of the number of transferred electrons 
n can be studied with the DO formalism. To do this, 
we unravel ps(t) in terms of this continuous projective 
measurement: p$(t) = X^n^s^W? similarly to how this 
is done in quantum optics [17]. The probability distri- 
bution of having n transfers after time t is given by 
P(n,i) = Trs{pg n ^(t)}, and the corresponding CGF is 

F(x\t) = m Er=-oc p (M) einx - This allows to calcu " 

late the A:-th order cumulant of the current distribution as 
(I k (t))c = {n k (t)) c = f t ^M^\ x ^ . In practice, count- 
ing in lead a can be effected by adding Xa to the tun- 
neling Liouvillian Ct through the replacement Vkam ~^ 
V kam e ipX( * /2 [18 , where p = +/- is the Keldysh in- 
dex corresponding to the forward/backward time branch. 
Derivatives with respect to different counting fields, e.g. 
Xio Xr? allow us to obtain also cross correlations of cur- 
rents flowing through different contacts. In the follow- 
ing, the lead-dependence of the counting field will be 
considered implicit. Let us try to relate this CGF (or 
alternatively the moment generating function Q = e^) 
with a general NM evolution. In the x-space, the density 
operator p s (x,z) = E5T=-oo Jo°° dtp^ (t)e ixn ~ zt follows 
the evolution ps(x, z ) = ^(x?^)Ps(0), with Q(x,z) = 
[z-W(X:z)}-\ and W(x, z) = C s + £(*,*). To lowest 
order we have 

S( X ,z) = Tr R |£ V ( X ) — J— ^- r v ( x )^| . (7) 
For later use, we also introduce the two-point self-energy 
n(x2, Xi, *) = Tr R |>Cv(X2) z _ c l_ £R ^v(Xi)PR q | -(8) 

Obviously, we have n(x, X: z ) = £(x, 2), and E(x = 
0, z) = E(z). Explicit expressions for Eqs. ^ and ([8| 
are given in appendix]] 

In the widely used Born-Markov approximation, the 
state at which counting begins (say t = 0) can be taken 
to be ps(0) p^. However, to consider NM corrections, 
the state at time t = can no longer be considered as 
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a separable state, as it contains initial system-bath cor- 
relations. To account for these, we explicitly divide the 
time evolution into two intervals (see Fig. [T]). The evo- 
lution from to (time at which system and reservoirs are 
separable) to t = is given by Z J_ C , while the evolution 
from t = to t is given by z _c( x ) • Doing this we obtain 
the moment generating function (MGF): 



(9) 



Here z is the conjugate frequency to t, and zq to —t^. We 
will take to —> — oo, which implies zo —> 0~ (henceforth 
implicit). The trace in ([9| refers to the full trace (system 
plus bath degrees of freedom). Using geometric expan- 
. of - — and and performing the trace over 



-C(x) z -C 

the reservoirs, we get 



l 



z-Cs- T,(x,z) 



(l+T( X ,z)) 



(10) 



In this equation, (...) = Trs {. . . p| tat }, where we have 
taken ps(0) = p| tat , as we are interested in fluctuations 
around the stationary state. This can be obtained ei- 
ther as p| tat = \im z ^o zps(z) m equation Q, or solv- 

AA</^\ ofof ^ -T 

mg 

Eq. (|10|) is given by 



W(0)p| tat = 0. The inhomogeneous term T(x,z) in 



T( X ;z) = i{n(x,o,z )-n( X) o^)} 



Eqs. ( 10 ) and ( 11 ) are the first main formal result of the 



paper. As we shall show below, the inclusion of T(x, z) in 
the MGF is crucial to account for NM physics and quan- 
tum noise. Importantly, T(x, z) cannot, in general, be 
cast in the form of a self-energy, since only one of the two 
vertices (i.e. tunneling Liouvillians) contains a count- 
ing field X- Notice that Eq. ( pTj ) extends the particular 
form of the inhomogeneity r(x; z) = |{E(0,0)-£(0,z)}, 
which appears in [19 j. This is only valid for a system 
with NM dynamics but with Markovian coupling with 
the bath in which counting is performed, and as a result, 
quantum fluctuations due to the Fermi contacts are not 
captured in this case. 



Noise spectrum 



From the MGF (10), together with (11), we can derive 
a general equation for the noise spectrum. To this end 
we make use of the MacDonald's formula [TOJ [20] 

/>OC 

5 (2) (w) = u / dtsm(ajt)(I 2 (t)) c 
Jo 



u 2 d 2 

' 2 d(i X y 



[Q(X, z = iu>) + (w ->• -w)] 



(12) 



(11) and obtain 



S {2) (u) = [(J n (iuj,iuj )) + { J\iuj, 1^)0,0(100) J 1 (iuo.iuoo))] + (u <-> -u), 



with ojq — >• and 



J H (z,z ) 
J\z,z') 



d 2 



d(ix2)d(ixi 
d 



-^(X2,xi,z) 



d(ix) 



X2,Xi^0 

[IL(0,x,z)+Il(x,0,z')} 



X-s-0 



r 



(13) 

(14) 
(15) 



Eq. (13), together with (14) and (15), is the second main 



formal result of the paper. It is exact and agrees with 
previous approaches in the literature in the appropriate 
limits [HJ [21]. In particular, the Markovian result [TO] 
is recovered by neglecting the frequency dependence 
of the jump super-operators: J JI (z, zq) —> J7 /J (0,0), 
J\z,z') J J (0,0). The correct NM zero-frequency 
limit [TO] is also recovered. It is interesting to notice 



that Eq. (10) not only allows us to obtain the NM noise 



spectrum, but also single-time NM correlations to arbi- 
trary order, (I N (t)) c , (n N (t)) c , (Q N (t)) c , by simply tak- 
ing derivatives with respect to the counting field. 



We notice that the above derivation has focused on 
particle currents flowing through the barriers separating 
central system and leads. At finite frequencies this parti- 
cle current is not conserved due to charge accumulations 
in the system, and the total current (particle plus dis- 
placement) needs to be considered to obtain the noise 
spectrum. However, our results are general, and current 
conservation can be considered by the inclusion of the 
proper counting fields [TO] in z) and Ti(x2->Xi-> z ). 
Thus, particle, total, and charge noise (equivalently volt- 
age noise for a capacitive system), can be calculated from 



Eq. ( 13 ). To this end, it is enough to consider respectively 
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E2] Xl/XR, Xtot = XL + Xflj and Xaccum = PXL ~ &XR, 

giving rise to different jump super-operators [TO] . 



RESULTS 
Single resonant level model 

We now use the formalism presented in the previous 
section to calculate the NM noise spectrum of a single 
resonant level model (equivalently of a single electron 
transistor with Ec ^> fceT, being Ec the charging en- 
ergy, and with only two relevant charge states) . The noise 
and charge spectrum of this system have already been 
calculated with a variety of techniques [3j [TTJ [23], and 
the exact solution is also well known [24]. We therefore 
use this as a benchmark of our method. In the follow- 
ing we show the good agreement between our theory and 
the exact solution. In the QNR, these two, in contrast 
to the Markovian result, show quantum-noise steps due 
to vacuum fluctuations, as we will see. The Markovian 
and non- Markovian results we present here correspond to 
first order in perturbation theory (sequential tunneling) 
and in the following S^ 2 \uj) refers to the 'total' noise. 

The SRL model is described by the Hamiltonian 



n = 



k,a£L,R 



E ^ c LlO)(l|+H.c. 



(16) 



k,aeL,R 



Here, each of the terms corresponds to central system, 
reservoirs, and tunneling respectively. The state |1) (oc- 
cupied level), together with |0) (empty level) form the 
Hubert space of the central system (|0)(0| + |1)(1| = 1). 
This model, despite its simplicity, contains a great deal of 
interesting physics: In the context of mesoscopic systems, 
this Hamiltonian captures the physics of a quantum dot 
in which only one single level participates in transport 
(strong Coulomb Blockade regime). Also, it can be shown 
that there is an exact mapping between the SRL model 
and the spin-boson model (namely a quantum two-level 
system coupled with strength a to an Ohmic dissipative 
bosonic bath) at a = 1/2. This mapping is actually 
an special case of the more general relation between the 
spin-boson model and the anisotropic Kondo model, for 
which a = 1/2 is the exactly solvable point, the so called 
Toulouse limit of the Kondo problem [25] . 

Fig. [2^1 shows the shot noise spectrum S^ 2 \uj) of the 
total current through the system obtained with the non- 
Markovian formalism discussed in the previous section 
(blue dashed-dotted curve). We also plot the exact re- 
sult [24] (black dotted curve) and the one obtained af- 
ter a Markovian approximation [10] (red dashed curve). 
The agreement between the exact solution and the NM 



calculation is extremely good. Both develop dips at fre- 
quencies uj = =b|e =b ^|, and show a strong frequency 
dependence. As expected, and due to the mapping afore- 
mentioned, the shot noise spectrum in Fig. [2^i agrees well 
with the one of a non-equilibrium Kondo model in the 
Toulouse limit [26]. In stark contrast, the Markovian so- 
lution is markedly different: it is frequency-independent 
and equals S^ 2 \uj — » 00) = 2 (rv!-f R ) ~ ^2~ m Even at 
uj = 0, the MA deviates from the NM and exact solu- 
tions, which here fall practically on top of each other. 
In Fig. [2}d, we explore the linear-response regime when 
the level is outside the bias voltage window. In this sit- 
uation shot-noise is negligible, and quantum fluctuations 
are dominant in the spectrum for Huj ^> fc#T. The quan- 
tum noise step expected at uj = e is fully captured by our 
NM approach, while here it becomes clear that the MA 
does not capture quantum noise physics. 

The richness of the SRL model can be further explored 
by noting that it also describes the physics of a single elec- 
tron transistor (SET) with charging energy Ec ^> &#T, 
and voltage such that only two charge states \N) and 
\N + 1) are relevant. One can describe a SET in this 
regime with Eq. ( fl6| ) by just making the substitutions 
[27] e E c , |0) \N) and 1) |JV + 1). Let us derive 
the charge- noise spectrum (2| of the SET. This prob- 
lem has already been studied by Johansson et al. using a 

different formalism [23] . As discussed in the previous sec- 

(2) 

tion, Sq (uj) can be found by considering the jump oper- 
ators arising form the counting field Xaccum = Pxl—^XRi 
being a and (3 coefficients determining how the total cur- 
rent is partitioned between both left and right contacts 
[TO] . Alternatively, we can apply charge conservation: 
Jl(£) —In(t) — Q(t), being Il/r the current through the 
left/right lead and Q the charge inside the well. This, 
together with the Ramo-Shockley partitioning theorem 
I(t) = aI L (t) + /3I R (t) to obtain 

1 



S { q\") = - 2 \St\") + ^» " " SrIH -(17) 

The cross correlations 

?(2) 



?(2) 



( 2 )/ 



?(2)/ 



s 



LR/RL 



/oo 
dre- <WT ({/ L/fl (r)/ B/L (0)}> c , (18) 
-CO 



can be easily calculated taking the derivative of the CGF 
with respect to counting fields xl and xr> while the 

(2) 

particle-noise contributions S L j R involve a double deriva- 
tive with respect to Xl/xr °f the CGF. Fig. shows 
the noise associated with the charge fluctuations in the 

(2) 

central island of an SET, Sq (uj). Interestingly, if the 
SET is used as a detector of another quantum system, 
this noise governs the measurement backaction p~3j [27] . 
When huj > Ec, the charge- noise spectrum contains 
extra quantum noise contributing to backaction, in full 
agreement with previous calculations [27] [28] . 

In Fig. [3] we investigate this zero-frequency limit given 
by our NM theory. Fig. [3^i shows the particle noise 
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FIG. 2: (color online). Quantum noise spectra of the SRL model (Tl = = 1 in all figures), a) S^ 2 \uj) as a function of 
frequency uj in the shot noise regime (e = 20, V = 100, T = 4). In this limit, the noise develops dips at uj = ±\e =b ^| . 
b) S^ 2 \uj) as a function of frequency uj in the quantum noise regime (e = 10, V — > 0, T = 1). In this limit, S^ 2 \uj) develops 
a quantum noise step at uj = e. c) Charge noise Sq\uj) as a function of frequency uj of a single electron transistor acting as a 
detector (Ec — V— 10). When uj > Ec, Sq\uj) contains extra quantum noise contributing to backaction. 
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FIG. 3: (color online). Zero-frequency limit of the non- 
Markovian theory, a) Particle-current noise for Tl — Tr, 
s/Tr — 80, T/Tr = 4 The noise suddenly increases when the 
level enters the voltage bias window, b) (0) as a function 
of bias voltage V for Tl =Tr, s = 0, T/Tr = 4. While the 
Markovian approximation is flat at all voltages, the NM and 
exact solution show a structure that strongly differs in both 
for low voltages due to cotunneling processes. 



S R (uj = 0) as a function of voltage for a configuration 
such that e/Tn = 80 ^> T/Tr = 4. We observe a res- 
onant step in the noise spectrum at precisely V = ±2s. 
Above this step, there is a discrepancy of the Markovian 
solution with the NM and exact results, while right be- 
low the step, Markovian and non-Mar kovian limits differ 
from the exact solution. This last discrepancy is due 
to cotunneling contributions, only captured by the ex- 



act result. The difference is better observed in Fig. [3}}, 
where we set e = and vary the bias voltage again. Re- 
markably, the Markovian solution is flat for all voltages, 
while both NM and exact solutions show certain struc- 
ture capturing system-bath memory effects. Only for low 
voltages these two disagree, when cotunneling contribu- 
tions become important. At zero voltage, the Marko- 
vian and NM curves coincide as expected (since the only 
contribution to noise should originate from equilibrium 
fluctuations). For large enough voltages, the exact and 
NM results fall on top of each other, and we remark that 
the limit V —> oo is exact in both Markovian and non- 
Mar kovian approaches, and thus all three curves converge 
to the same value in this limit. 



Single-time full counting statistics 



Beyond frequency-dependent noise spectra, Eq. (10) 



also allows us to study single-time full counting statistics 
of the number of electrons n transferred to a particular 
terminal. This quantity is denned through the cumulant 
generating function T = log Q as 



(n k (t)) c 



gk 



d(i X Y 



Hx\t)\ 



x=o- 



(19) 



Such /cth-order cumulants can be measured by e. g. 
counting electrons using a quantum point contact and 
analyzing the time-dependent statistics of the events 
[29]. Fig. [i^i shows the single-time Fano-factor = 
( n i?W)c/ ( n R) °f the SRL model. This figure shows up 
to the fifth order (k = 5) Fano-factor (solid lines) to- 
gether with the results corresponding to the MA (dotted 
lines) in the shot noise regime (level within the bias volt- 
age window). At large times, the agreement with the 
Markovian solution is good for all the Fano-factors. At 
short times, however, the MA converges to the Poissonian 



FIG. 4: (color online). Single-time full counting statistics 
of the SRL model (Fano-factors of the right particle-current 
F (k) (t) = (n k R (t))c/{n R ) with F L = F R = 0.25, T = 1) up to 
the fifth cumulant. a) Shot noise regime (e = 20, V = 100). 
b) Quantum fluctuations regime (e = 20, V = 30). Dotted 
lines correspond to the Markovian approximation. 

limit while the NM solutions clearly show a strong sub- 
Poissonian supression. More interesting is when the level 
is above the bias window and all noise comes from quan- 
tum fluctuations (Fig. ^p). In this case, and taking an 
infinite bandwidth, the second cumulant C2(t) = (n 2 (t)) c 
can be approximated as the inverse Laplace transform of 

with \I/ the digamma function. This gives the expo- 
nentially large Fano-factor ^(t) « for very short 
times, and follows the power law 1*2 (£) ~ ^{e-^id 
at intermediate times. From this result we can estimate 
the time at which F^(fy deviates from the MA, namely 

Switch = ^ £ -^ R )i dc ' In Fi §-0 3 we P lot this P° wer law be- 
havior (dashed blue line) together with the full NM solu- 
tion (solid lines), and the Markovian solution, which here 
lie at the Poissonian value 1. For times t <C Switch? we 
obtain large super-Poissonian noise resulting from high- 
frequency quantum fluctuations. 

Double quantum dot 

To further illustrate the theory, we now consider the 
example of a double quantum dot (DQD). To the best 
of our knowledge, a complete study of this model in 
the different regimes of V, T and cj, and in the NM 
limit is yet lacking. The following results are also ap- 
plicable to a Cooper pair box qubit. Again, the Marko- 
vian and NM solutions shown here correspond to first 
order in perturbation theory (sequential tunneling) and 
S^ 2 \uj) refers to the 'total' noise. In the Coulomb block- 
ade regime, the possible DQD states are |0) = \Nl,Nr), 
\L) = \N L + 1,N R ) and \R) = \N L ,N R + 1), with N L /N R 
being the number of electrons in the left /right dot. The 
qubit, with Hamiltonian Hs = e (\L)(L\ - \R)(R\) + 
T C (\L)(R\ + \R)(L\), has eigenvalues E± = ±f, being 



FIG. 5: (color online). Quantum noise processes in a dou- 
ble quantum dot. In the QNR, quantum fluctuations can 
discharge the system through the left/right reservoir if ftw > 
\/jll/r — A/2\. These correspond to the steps in Fig. [6^l. When 
uj — A, quantum interference between the eigenstates gives a 
noise suppression. 



A = 2^/ e 2 + T c 2 . Near linear response (eV <C &bT, Huj), 
the only noise contribution originates from equilibrium 
fluctuations - either thermal noise for ksT ^> hw, or 
quantum noise for huj ^> ksT. In Fig. [5] we sketch the 
physical processes due to quantum fluctuations, which 
give rise to the noise spectrum in Fig. [6^i. For eV < A, 
the conductance is zero and therefore S^ 2 \0) = 0, as 
dictated by the fluctuation-dissipation theorem. Quan- 
tum fluctuations, on the other hand, give rise to a finite 
noise for uj > (steps at Huj = \[1l/r ± ^| in Fig. 6a). 
Importantly, this physics is not captured with the MA, 
neither by other models for the inhomogeneity, such as 
T(X5 Z ) = iH^(XjO) — E(x, z)}. The spectrum also con- 
tains a strong dip centered at uj = A. This dip, which is 
voltage-independent and reaches S^ 2 \uj = A) = 0, can 
be understood as resulting from coherent destructive in- 
terference between the qubit eigenstates. This is demon- 
strated in Fig. (6}d, where we investigate how this feature 
at uj = A changes as we move the Fermi energy, Ep, of 
the reservoirs. For V = 0.1 and Ep = (black solid 
curve), E+i is above/below the chemical potentials and 
we find a dip shape, as discussed. When Ep is aligned 
with the lowest level, namely Ep = E = — ^, the res- 
onance changes to a Fano shape, as one expects from 
interference between a discrete level (the one above the 
chemical potentials at E+ = ^) and one strongly cou- 
pled to a continuum (the one at Ep = E = — y). When 
both levels are above E Fl the interference at uj = A is 
suppressed (red dotted curve). However, if both levels lie 
above Ep (light grey curve), quantum interference still 
occurs, giving in this case a narrow resonant peak in the 
noise spectrum, since now we have a qubit weakly cou- 
pled to the leads - therefore with a low dephasing rate. 
A very important remark of this figure, is that the situ- 
ation corresponding to Ep = — 4 gives a different result 
from that corresponding to Ep = 4. In the former, the 
peak at uj = A has been suppressed, while in the last, 
the resonance occurs. This we understand in terms of 
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coherent oscillations only taking place when the levels 
lie below the chemical potentials. Most importantly, the 
light-grey curve only presents one quantum noise step, 
corresponding to the anti-bonding state. As the charge 
oscillates fast between both eigenstates, this can decay 
to the reservoirs via quantum noise processes only from 
the lowest level. However, in the situation with both 
eigenstates above the chemical potentials, charge can de- 
cay to the reservoirs from both levels through quantum 
noise processes. If eV > A, transport is possible and 
shot noise is finite, therefore *S^ 2 ^(0) 7^ 0. This limit is 
discussed in Fig. [6)3. Interestingly, quantum noise is pro- 
gressively overcome by shot noise as V increases. As a 
result, for large voltages, the quantum noise steps dis- 
appear and the noise is of smaller magnitude. In this 
case an incomplete destructive interference is found at 
uo = A: S^ 2 \uj = A)/Idc(V) is greater than zero and 
does not depend on V. The width, on the other hand, 
increases with the voltage, which can be understood as 
a decrease of the dephasing time (inverse of the width) 
due to the coupling with the reservoirs [5]. The MA is 
recovered as V — >• 00, with features at uo = and uo = A 
on top of a background of sub-Possonian partition noise, 
Fano-factor (uo)/I dc (V) = 1/2. 

The transition from a Fano shape to an anti-resonance 
in the noise spectrum encountered in Fig. [6}d is fur- 
ther investigated in figure [7^i. Here we show how the 
quantum noise step progressively appears as the bonding 
state comes below the chemical potentials. At the same 
time, the Fano resonance gives rise to the destructive- 
interference feature at the qubit frequency. The effect of 
temperature is shown in Fig. [7}}. Still in the linear re- 
sponse regime, where the 'shot' contribution is negligible, 
we see how quantum noise is overcome by thermal noise, 
giving a finite value at zero frequency for increasing 
temperature, as dictated by the fluctuation-dissipation 
theorem. The Fano shape, consequence of having the 
lowest level strongly coupled to the reservoirs, but also 
coupled to the anti-bonding state, persists at high tem- 
peratures. 

CONCLUSIONS 

We have presented a general non-Markovian theory of 
frequency-dependent noise based on QME. The impor- 
tance of NM correlations to correctly capture the physics 
of vacuum fluctuations has been shown through the study 
of a single resonant level model and a double quantum 
dot in the quantum noise regime. Our equations for the 
CGF and noise spectrum open the possibility to investi- 
gate this physics in a variety of systems where NM correc- 
tions are of vital importance, such as electromechanical 
resonators close to the zero-point motion [30] , or strongly 
correlated cold atoms in optical lattices [3T] , 
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Liouvillian Perturbation Theory 

For concrete application of the formalism, we employ 
Liouvillian perturbation theory (LPT), as described in 
Refs. [T4HT6] . We give here a brief review of the essential 
elements of this theory — for more details, the reader 
is referred to the original references. As explained in 
the main text, the Hamiltonian describing single-electron 
tunnelling between central system and reservoirs is 

U T =J2 V kamcl a d m + V£ am S m C koi , (21) 
kcxm 

where d m is the annihilation operator for the single- 
particle level m in the central system, and c ka the an- 
nihilation operator for an electron with momentum k 
in lead a. and Vkcxm is a tunnelling amplitude. Intro- 
ducing a compact single index "1" to denote the indices 
(£1, fci, ai), we have 

ci = c 6lkiai = \ fZl > ( 22 ) 

L c k 1 a 1 •> Si — 

with index £i = ± indicating whether the operator is a 
creation or annihilation operator. We further define the 
system operators g ka = J2 m Vkamdm, such that the tun- 
nel Hamiltonian can be simply written as Ht = 
where 1 denotes (— £i, fci, ai), and here as elsewhere, im- 
plicit sums over repeated indices. In the same notation, 
the reservoir Hamiltonian H res = ^Zk,aeL,R e kac\ (X c ka 
reads H res = £iCic T % + , 

In Liouville space the tunneling Liouvillian can be writ- 
ten as 

C T = -i[H T ,.} = Y, V° V C\G\, (23) 

P 

where p = ± is a Keldysh index corresponding to the two 
parts of the commutator. C and G are the superoperators 
corresponding to c and g, defined through their actions 
on arbitrary operator O: 

c p = r oio, P = + (24) 

1 [ Oai, p = - 

G\0 = a p x { P = + . (25) 

1 I -Ogi, p = - 

The object <j p is a dot-space superoperator with matrix 
elements (cr p ) ss , = 5 S s5 s 's'P Ns ~ Ns ' where N s is the 
number of electrons in state s. 
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FIG. 6: (color online). Finite-frequency noise of a double quantum dot (results normalized to the dc current in the large- 
voltage limit I dc (V = oo) = t 2 c F R /[F 2 R /4 + 4s 2 + t 2 c {2 + r^/r^)]. a) Near linear response, S^ 2 \uS) shows quantum noise steps 
at uj — |A/2 d= V/2| and a dip centered at uj — A (indicated with arrows for the case V — 4 in the figure), b) The feature at 
uj = A originates from quantum interference between the bonding and anti-bonding qubit states. Here we set V = 0.1 and vary 
the reservoir Fermi energy Ef, observing a displacement of the quantum noise step, as well as a modification of the resonance 
form at the qubit frequency A (see text), c) Shot noise limit. Quantum noise steps are only visible for V < A, fcsT, otherwise 
the contribution from shot noise or thermal noise are dominant. Parameters: e — 0, A = 2T C = 6, Tl — Fr = T/2 = 0.1. 



We can now express the self-energy in terms of these 
superoperators. As described in the main text, the non- 
Markovian system self-energy can be obtained by expan- 
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FIG. 7: (color online), a) Effect of a gate voltage. As the 
relative distance between the dot levels and the lead chemi- 
cal potentials is varied (here illustrated decreasing the Fermi 
energy Ef), a quantum noise step, absent when the bond- 
ing state is aligned with both chemical potentials, appears 
at the corresponding frequency difference. The Fano shape, 
however, gives an anti-resonance at the qubit frequency A. 
Here T = 0.2. b) Effect of the temperature. As T is in- 
creased, the quantum noise step is lost, since thermal noise 
overcomes quantum noise, giving a finite value at zero 
frequency. The Fano shape is however preserved for high tem- 
peratures. Here Ef = —2.5. In both figures V = 0.1, e = 0, 
A = 2T C = 6, F L = F R = 0.1. 



sion of 

ps(z) = Tv res {[z - C}- 1 pg(t ) ® = n(z)p s (t ), 

as a power series of Ct and tracing out bath degrees 
of freedom. This can be done using the diagrammatic 
technique explained in Refs. [T4HT6] . To lowest order 
(sequential), we obtain 



(26) 



In this expression we find the free-propagator fts(z) 



z—Cs 



, and the reservoir contraction 



121 



= %Pi/(-£iPi£i), 



(27) 



with Fermi function f(e a ) = (e £_/ia + 

Counting in lead a is introduced through the replace- 
ment Vkam — ^ Vkame %PXa ^ 2 in the tunnel Liouvillian Ct- 
The x-dependent self-energy is then simply obtained as 
the above self-energy but with x-dependent superopera- 
tor G(x) replacing G. The two-point self-energy deter- 
mining T(x, z) can similarly be derived. We obtain 

nOfe.Xi,*) = -27ipiGf \4> a )){{cj> a \G^ 

Xe *^«l^l(PlXl-P2X2) 

xlgiAa + b^-iz), (28) 

with A a , 1 a )), eigenvalues and eigenvectors of the 
central-system Liouvillian, that is Cs\^ a }) = — zA a |0 a )), 
and Ip 2 \\) defined as 

/M 



did 



z0+ + puj — A ' 



(29) 



The upper limit of this integral can be taken as a 
Lorentzian cutoff D(uo) = X^/(uj 2 + X%), which gives 



J?>W = ^){i/W + ^)-^}. 



(30) 
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Here, 0(A) = \ ( g(X) + g(-\) - 2g(iX c )) and g(u>) = 
"3/ (| + ^i), being ^(x) the digamma function. In the 
wide-band limit (X c ;§> to, A), this integral becomes 



1 



ip 



/(pA) + ^(A), 



with approximate ^-function 

0(A) = J (5(A) - g(-X)) - log(X c /27r). 



(31) 



(32) 



This latter result is adequate for finite-frequency shot- 
noise calculations, but the more accurate form Eq. (30) 



is required to correctly capture the the single-time full 
counting statistics, for which a bandwidth X c — 500/c#T 
was assumed. 
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